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Abstract: We study scalar solitons on the fuzzy sphere at arbitrary radius and noncom- 
mutativity. We prove that no solitons exist if the radius is below a certain value. Solitons 
do exist for radii above a critical value which depends on the noncommutativity parameter. 
We construct a family of soliton solutions which are stable and which converge to solitons 
on the Moyal plane in an appropriate limit. These solutions are rotationally symmetric 
about an axis and have no allowed deformations. Solitons that describe multiple lumps on 
the fuzzy sphere can also be constructed but they are not stable. 
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1. Introduction 

It is an interesting fact that the classical equations of motion of noncommutative field the- 
ory have a richer set of solutions than their commutative counterparts. Noncommutative 
solitons include deformations of conventional solitons but also new objects that are absent 
in the commutative theory. This applies to both gauge theories and scalar field theories 
and can be traced to the existence of a new length scale associated with the noncommu- 
tativity. A striking example is provided by scalar theory on the Moyal plane [1], where 
the corresponding commutative theory can have no solitons by Derrick's theorem [2]. Sev- 
eral authors have found explicit soliton solutions in gauge theory with and without matter 
fields, see, e.g., [3-6] and noncommutative deformations of commutative solitons have been 
studied in [7-9]. For reviews of noncommutative quantum field theory we refer to [10,11]. 

The scalar solitons on the Moyal plane were first observed in the limit of infinite 
noncommutativity where the kinetic term in the action can be neglected and solutions 
correspond to multiples of projectors [1]. This line of reasoning was generalised to scalar 
fields on arbitrary Kahler manifolds including the sphere in [12] and subsequently multi- 
solitons on the fuzzy sphere in the infinite noncommutativity limit were considered in [13] . 
In [14,15] it was shown that there exist scalar solitons on the Moyal plane at large but finite 
values of the noncommutativity parameter and that no soliton solutions exist depending 
smoothly on small values of the noncommutativity parameter as expected from Derrick's 
theorem. All the solitons at finite noncommutativity constructed so far have the property 
that they are rotationally invariant about some point. It is not known whether there 
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exist stable soliton solutions in the plane describing separated lumps, but there are good 
reasons to believe that no such solutions exist and that all soliton solutions in the plane 
are rotationally invariant about some point. First of all, separated lumps attract each 
other according to perturbative calculations at large noncommutativity parameter [16,17]. 
Furthermore, it has recently been proven that there cannot exist a family of static solutions 
on the Moyal plane at finite noncommutativity that interpolates smoothly between the 
solution describing two overlapping solitons and a solution with two infinitely separated 
solitons [18]. Finally, the stability results for scalar solitons on the fuzzy sphere obtained 
below, combined with a scaling limit that yields the Moyal plane from the fuzzy sphere 
(see e.g. [19]), provide further evidence for the non-existence of stable multi-lump solutions. 
Further results for noncommutative theory in the plane can be found in [20-25], and on 
the fuzzy sphere [26-33]. 

In the present paper, we study the equations of motion for a scalar field on the fuzzy 
sphere in detail, paying particular attention to issues that arise at finite noncommuta- 
tivity. Part of our analysis is parallel to that of [14, 15] which dealt with existence and 
stability of scalar solitons on the Moyal plane. The algebra underlying the fuzzy sphere 
is finite dimensional, and this simplifies the mathematics, especially the stability analysis. 
It enables us to establish some uniqueness results by ruling out deformations of the field 
profile of a stable soliton. We are also able to prove a stronger nonexistence result at small 
noncommutativity than is available on the Moyal plane. 

The fuzzy sphere is characterised by two parameters: R which can be interpreted as 
the radius of the sphere, and j which labels the associated matrix algebra. Together these 
parameters determine the strength of noncommutativity. A scalar field theory on the fuzzy 
sphere is further characterised by a scalar field potential V(<p) which we take to be a double 
well potential with a global minimum at the origin. Our main results are then the following: 

• There are no solitons for small values of the radius. 

• As the radius is increased a large number of soliton solutions come into existence, 
most of which are unstable. 

• A necessary condition for a soliton, that is rotationally invariant around the north- 
south axis, to be stable is that its eigenvalues form a decreasing sequence (in the 
standard su{2) basis). At sufficiently large radius, another condition is that the 
eigenvalues lie close to the minima of the potential, and these two conditions are 
also sufficient. 

• The only small deformations of a stable soliton which lead to new soliton solutions 
are rotations. 

• There is a family of stable, rotationally invariant solitons on the fuzzy sphere which 
converge to stable solitons in the Moyal plane as both R and j tend to infinity keeping 
R 2 /j fixed and sufficiently large. 

The limit of infinite noncommutativity is by no means smooth, and the set of classical 
solutions is drastically reduced at finite as compared to infinite noncommutativity. The 
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difference is particularly striking when one considers multisolitons. In fact, the very notion 
of a multisoliton solution is somewhat problematic at finite noncommutativity since there 
is no precise definition of soliton number in this case. In the limit of infinite noncommu- 
tativity, the soliton number of a stable classical solution can be taken to be the rank of 
the corresponding projector, but solutions at finite noncommutativity are no longer given 
by projectors and this definition fails. A physicist might wish to define a multisoliton as 
a configuration with two or more localised lumps in the scalar field that are separated on 
the fuzzy sphere, but this is not a precise mathematical definition. Our results strongly 
suggest that there are in fact no stable solutions with separated lumps, and that all sta- 
ble solutions are rotationally symmetric. There is, in other words, no moduli space of 
scalar multisolitons on the fuzzy sphere at finite noncommutativity. We do, however, give 
evidence for multi-soliton solutions on the fuzzy sphere describing separated lumps that 
are evenly distributed along a great circle. These configurations are unstable due to the 
mutual attraction between lumps. In the special case of two identical lumps at the north 
and south poles, we give rigorous arguments for the existence of the unstable solutions. 

By taking an appropriate limit j, R — > oo, one can blow up the region close to the south 
pole into the Moyal plane. Assuming that every solution in the Moyal plane descends from 
solutions on the fuzzy sphere in the limit, one would expect a solution with separated lumps 
on the Moyal plane to descend from an infinitesimal deformation that breaks the rotation 
symmetry of a solution on the fuzzy sphere. The fact that no such deformations exist 
strongly suggests that there are no stable multi-lump solutions at finite theta on the Moyal 
plane, in agreement with the results of [18] and expectations from perturbative calculations 
in the limit of infinite non-commutativity [16,17]. 

The unstable great circle solitons are lost in the Moyal plane limit since any lump 
not sitting at the south pole will be moved to infinity. By symmetry we would, however, 
expect the existence of an unstable solution consisting of an infinite chain of lumps finitely 
separated in a straight line on the Moyal plane. Such a solution would be periodic in one 
direction and have infinite energy. It could presumably be further generalized to a regular 
two-dimensional array of lumps. 

The paper is organized as follows. In the following section we give a brief introduction 
to the fuzzy sphere and establish our notation. In section || we prove the nonexistence 
of static soliton solutions at small radius. In section |] we study unstable solitons which 
bifurcate from the trivial soliton as the radius is increased. We construct some rotationally 
symmetric solitons on the fuzzy sphere in section || show that they are stable in section |6|, 
and that they converge to solitons on the Moyal plane in an appropriate limit in section |8|. 
As a byproduct of the stability analysis in section ||, we show that the only possible defor- 
mations of these stable solitons correspond to rotations that do not affect the shape of the 
field profile. In section ^ we establish the existence of (unstable) solutions that correspond 
to separate lumps at the north and south poles and give arguments for their generalization 
to lumps evenly distributed along a great circle. 
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2. The Fuzzy Sphere 



One of the simplest interesting noncommutative spaces is the fuzzy sphere [34] which can be 
regarded as an approximation to the true sphere obtained by cutting off the function algebra 
on the sphere by requiring the total angular momentum quantum number to be smaller 
than a prescribed value j. This requires of course that the ordinary product of functions 
be modified to a noncommutative product. If x%, X2 and X3 are Cartesian coordinates in 
M 3 and the sphere is defined as 

x\ + xl + xl = R 2 , (2.1) 

then the commutation relations between the Xi's become 

[xj,x k ] = iae jk ex£. (2.2) 

Here a is a positive real parameter. For our purposes it is convenient to regard the fuzzy 
sphere as the algebra Mj of (2j + 1) x (2j + 1) matrices, where j is a half integer. These 
matrices act on a (2j + l)-dimensional vector space which we denote by Wj. In Mj we 
have the spin matrices (su(2) generators) for spin j which satisfy the usual commutation 
relation 

[Jj, Jk] = if-jkiJi, (2.3) 

and JiJi = j(j + 1). These matrices can be seen to generate the full matrix algebra Mj. 
With the above conventions we see that the X{ are proportional to the angular momentum 
generators, Xi = aJi, and the noncommutativity parameter a is governed by the spin j 
and and the radius of the sphere through 

R 2 =j(j + l)a 2 . (2.4) 

The commutative sphere of radius R is recovered in the limit j — > 00, keeping R fixed. 
Alternatively, one can take a scaling limit j,R — ► 00, keeping the ratio R 2 /j fixed, and 
consider a small neighborhood around the south pole. The fuzzy sphere commutation 
relations fl2.2|) can then be shown to reduce to those of the Moyal plane [19]. 

Scalar field theory on the fuzzy sphere [35] is defined by the action functional 

S=——Tr {[Ji, </>][$, Jil+I&VM) (2.5) 

where V is the potential and 4>, the field, is an arbitrary hermitian (2j + 1) x (2j + 1) 
matrix. The action is invariant under rotations <f> — ► TZcpTZ^ 1 where 1Z = ex.p(i9hiJi). The 
variational equation of S is 

2[Ji,[Ji,cf>}}+R 2 V'((f>)=Q. (2.6) 
It follows immediately that any solution to this equation satisfies the condition 

TrV'(<f>) = 0. (2.7) 
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Introducing the raising and lowering operators J± = J\ ±iJ2, the variational equation can 
be written 

[J_, [J+,0]] + [J 3 ,0] + [J 3j [J 3 ,0]] + ^'(0) = 0. (2.8) 

Let \m) denote the standard basis for Wj consisting of the eigenvectors of J3, J^\m) = 
m\m), m = —j, — j + 1, . . . ,j. If a solution (ft to Eq. ( |2.g| ) is diagonal with respect to this 
basis, its matrix elements <f) m = (m\(j)\m) satisfy the equation 

R 2 

a m (4>m+l - 4>m) ~ a m -l(4>m, ~ 4>m-l) = -^^'(^m) , (2.9) 

for m = —j, — j + 1, . . . , j, where a m = j(j + 1) — m(m + 1). It is convenient to formally 
introduce <f>j+i and <ft-j-i with the convention that (fij+i = <A? an d = tft-j- Summing 

eq. ( |2.9| ) over m the left hand side telescopes and we obtain the first order difference 
equation 

R 2 m 

4> m+ i - 4> m = ^2v'(4>i), m = -j,---,j-l (2.10) 

l=-J 

with the constraint 

E V'^i) = 0. (2.11) 
i=-j 

Diagonal matrices commute with J 3 and diagonal solutions of the scalar field equations 
correspond precisely to solitons which are invariant under rotations about the z-axis. 

We shall restrict ourselves to studying potentials V(x) which are twice continuously 
differentiable, nonnegative, and having a double zero at x = and V(x) > for i^O. We 
shall assume V(x) has a local minimum (false vacuum) at s > in addition to the global 
minimum (true vacuum) at x = 0, that is V'(x) has zeros at x = 0, x = r and x = s with 
< r < s, and no others (see figure [l]). 

We first note that any solution <fr to the equations of motion ( |2.6D has its spectrum in 
the interval [0, s]. A proof of the equivalent result for scalar field theory on the Moyal plane 
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is given by the authors of [15]. This proof carries over immediately to the fuzzy sphere so 
we do not reproduce it here. 

If <j> is any operator on Wj we denote its matrix elements (m|^>|n) by <p m ,n- Before 
proceeding, we note that, in addition to rotational symmetry, the soliton equation ( p.6[ ) is 
invariant under the transformation 

^m,n -> <t>-m-n ( 2 -12) 

which corresponds to reflection in the equator. Then, for example, if is a diagonal soliton, 
so is ifi, where tp m = 4>- m . 



3. Nonexistence at small R 



In this section we show that, at fixed j, there are no solitons when the radius R is suf- 
ficiently small. For convenience we define /i = R 2 /2. We define an inner product and a 
corresponding norm on hermitian matrices by 

{4>^) = Tt{U), \\<P\\ 2 = Tt4> 2 . (3.1) 



Let (ft be a solution to the soliton equation ( |2.6| ) and let / denote the unit matrix. We 
begin by showing that 

\\<f>-rl\\ = 0(n). (3.2) 



The notation / = 0(ip n ) will mean ||/|| = 0{ 

In order to prove fl3.2|) we note that any hermitian matrix <j> can be split into two pieces 

4> = <t>\\ + 4>± (3-3) 

with proportional to the identity matrix, and (<^ii, <?!>_l) = (so that Tr <^j_ = 0). The 
Laplacian operator (total angular momentum operator) in ( [2.6D has eigenvalues k(k + 1) 
with k = 0, 1, • • • , 2j. The unique eigenstate with k = is the identity matrix. Substituting 
4> into fl2.6| ) gives 

[JMM\+*v'<<i>) = v (3-4) 

so that 

M||^)|| >2||0j.||. (3.5) 
Since any solution <p has its eigenvalues in the interval [0, s], we find that 

\\<t>±\\ < (3.6) 

where c is a constant. Next, we look at the trace condition 

TrV%\+<f> ± ) = 0. (3.7) 

Then 

TrV'ih) + 0(fx) = (3.8) 
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and we conclude that for small values of ji 



011 =x Q J + 0(/x), (3.9) 

where xq is a zero of V', and so 

<f> = x Q I + 0(ji). (3.10) 

Next, we diagonalise <f> and discover eigenvalues 4>m = x o + Since every eigenvalue 

lies in the range [0, s], the trace condition implies that there must be some eigenvalues 
smaller than r and some larger than r. So by taking ju small enough, we must have xq = r 
and this proves equation (|3.2| ). 
If we now write 

4> = rl + x, X = X\\+X±, (3-11) 
using the splitting (|3.3[) , and substitute into the soliton equation fl2,6| ), we obtain similarly 



n\\v'(4>)\\ >2||x±||. (3.12) 

On the other hand we have 

||y'(^)|| 2 = Try'(r/ + X ) 2 

= ^'(r) 2 ||x|| 2 + 0( X 3 ) (3.13) 
= V"(r) 2 (\\x±\\ 2 +\\x\\\\ 2 ) + 0(x 3 ). 

The trace condition gives 

TrV'((j ) ) = Tr(x± + X\\)V"(r) + OU 2 ) (3-14) 

and since Tr x± = 0, 



|X||||=0(X 4 ) (3.15) 



so 



||y'(^)|| 2 = y"(r) 2 || X± || 2 + 0( X 3 ). (3.16) 
Comparing this with ( |3,12| ) we find 

M 2 (^"(r) 2 ||x±|| 2 + 0(x 3 ))>4||x±|| 2 . (3.17) 

Noting that | \x±\ | 2 = | |xl | 2 — | |x|| 1 1 2 = I Ixl | 2 (1 + C(x 2 )) we obtain finally 

/ u 2 (y"(r) 2 + 0(x 3 ))>4. (3.18) 

If // is small enough, this is impossible. 

Recall that at fixed j the radius R determines the non-commutativity parameter a 
through ( [2.4p . We have thus established the non-existence of scalar solitons at sufficiently 
weak non-commutativity on the sphere. This is certainly in line with previous results on 
the Moyal plane [14] but this was not guaranteed a priori given that Derrick's theorem 
only holds on the commutative plane but not on the sphere. 
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4. Existence close to critical values of R 



In this section we establish the existence of solutions to the diagonal soliton equations 
Q2.10 ) and ( 2,11| ) at certain critical values of the radius R. When we perform our stability 



analysis in section [6], these solutions will all turn out to be unstable. They do, however, 
provide insight into the structure of the field equations and we include them before moving 
on to the construction in section |||, which yields stable solutions for sufficiently large R. 

We think of (fi-j as the first entry in the matrix and then use the first order difference 
equation ( |2.10| ) to find each 4> m for m > —j as a function of 0_j = x. In order to discover 
whether a particular value of x gives a soliton solution, we only need to check the constraint 



( 2.11 ). This is equivalent to looking for zeros of the function 

g,(x) = V'Mx)] (4.1) 
i=-j 

where fj, = R 2 /2. We know that there are (trivial) solutions to the soliton equations at 
x = 0, r,s so g^x) must have zeros at these points. We also know that when [i is small, 
<? M (x) has no other zeros. Also, if x < 0, the (j) m {x) get locked into the negative region of 
V and so g^x) < 0. Similarly, when x > s, g^x) > 0. 

Assuming that solitons do exist at some values of we can imagine increasing [i to 
the critical point [i c at which they come into existence. Imagine plotting a graph of g^ix) 
against x. As we move through the critical value of /i, g^(x) will brush and then cut 
through the x axis. This will happen at [i c ,Xc satisfying 

g^c) = g'^(x c ) = o, (4.2) 

and by the continuity of g^,(x), we will generically pick up a pair of new solitons. 

In order to find the critical points, one in general needs to solve a nonlinear equation. 
There is, however, a series of critical points with x c = r which we can determine analytically. 

Setting = x = r, we solve the soliton equations and find (fi m = r for all m, so 
that <?«(r) = 0. Differentiating Eq. (|4.1| ) with respect to x and putting x = r leads to the 
constraint 

<(r) = V"(r) £ #(r) = 0. (4.3) 
i=-j 

We can differentiate (|2.9| ) to obtain a difference equation for the <^(r), with initial value 
tf-jir) = 1 

a m W m+ l{r) ~ <P' m {r)) - « m -i(C(r) - ^(r)) = ^V" {r)<P' m {r). (4.4) 

Let us define a diagonal matrix u with ua = (fc(r) (note that, although we have not 
explicitly written it, the matrix u is a function of //, inherited from the definition of the 
4>i). As long as /i 0, we can combine ( [4.3|) with ( |4.4| ) to obtain an equivalent matrix 
equation 

[J i ,[J i ,u]}+ f xV"(r)u = 0. (4.5) 
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The Laplacian has eigenvalues k(k + 1), with k = 0, 1, 2, • • • , 2j, and there is precisely one 
eigenstate represented by a diagonal matrix e k for each eigenvalue. Therefore the critical 
values of fi at x c = r are 

-k(k + l) , , N 

Ate = TFjr, s ^ , fc = 1, 2, • • • , 2i, (4.6) 



and at the critical point we have 



V"(r) 



it(// c ) = ae fc , a / 0. (4.7) 



We note that the number of critical points increases with j but the lowest critical value of 
\i (and therefore of R) is independent of j. 

The picture that we have is the following. As we increase \i towards a critical point, 
the function g(x) flattens at x = r, and as we move through the critical point, we expect g 
to start to cut through the x-axis on either side of x = r so that we gain two new solitons. 
Alternatively, it is possible that this process will precisely correspond with another pair of 
solitons moving in towards the critical point leading to a cancellation so that the net effect 
is the loss of a pair of solitons. In any case, when g'(r) changes sign, we must either gain 
or lose two soliton solutions since g(x) is continuous. 

There is an additional possibility which we cannot yet discount. Even though g'(r) = 
at the critical value of /U, it is possible that it may not change sign as we move through the 
critical value. This can only be if 

d_ 

dji 



■«£(r)U=*» = 0. (4-8) 



If we assume (4.8) is true, then by a very similar argument to the previous, we obtain the 
matrix equation 

[J h [JiM^M + V"[r)u{ji c ) + n c V"(r)u{ii c ) = (4.9) 
where u denotes u differentiated with respect to \x. As before, we can expand u(fi) in terms 
of the eigenstates of the Laplacian u(fi) = YlfLo u i(^) e i-- We- are studying a critical point 
satisfying [i c V"{r) = —k(k + 1) for some fixed k, so let us pick out the coefficient of e&: 

= k(k + l)«fe(A*c) + V"{r)u k {ii c ) + v c V"(r)u k (vc) = V"(r)u k (^ c ). (4.10) 

But by (|4.7|) , u{jji c ) = Uk{^ c )ek 7^ 0, so this is a contradiction. 

The conclusion is that whenever \x moves through one of the critical points (|4.6| ), we 
either gain or lose two solitons. As a corollary, there is a range of [i close to each of the 
critical points for which solitons exist. 

To illustrate, figure |2| shows an example plot of the zeros of g^{x) for j = | and a 
potential with V'(u) = u(u — 1.2) (u — 2.0). We see that there are many critical points in 
addition to those we have calculated. In most cases, soliton solutions are created as we 
increase fi, but they are also occasionally lost. In some cases, the critical points occur as 
bifurcations of existing solitons, but in most cases a pair of solitons appears (or disappears) 
where there was none before. Nevertheless, we call this the bifurcation diagram. We have 
observed that if one increases j, the complexity of the bifurcation diagram increases, and 
the number of solitons at sufficiently large fi increases. 
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X 



Figure 2: Diagram showing the zeros of g^x) for j 



5. Diagonal solitons 



We now show, for any given value of j, that the diagonal soliton equations ( 2.10| ) and ( 2.11| ) 
have a solution provided R is large enough. In section || we will further show that this 
construction leads to stable solutions. As j increases, we will need to increase R? linearly, 
so let us write R 2 /2 = 9j at the outset. When we come to consider the planar limit, it will 
also be convenient to rewrite the matrix elements of d> 



2/ 



^2 4>m\m)(m\ =^2^] - j +q)(-j +q\ 

q=0 



m=-j 



SO that ijjq 



U) 



-j+Q 



. Then the equations of motion ( 2.10| ), ( p. 11 ) become 



(5.1) 



2(g+l)(l-g/2j) 



k=0 



2j 



k=0 



(i) 

We have included an index j in ipq ' to remind ourselves that we are working with the fuzzy 
sphere, Mj, which is the j-th representation of su(2). 

We will see that solutions do exist for sufficiently large 9 by adapting the proof given 
in [14] for the Moyal plane. 



(5.2) 



(5.3) 



Theorem 1 

There is a constant, 9 C , independent of j such that when 9 > 9 C equations t\5.^j), ([5.^ have 
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a non-trivial solution. 



To see this, suppose V has its first local maximum in the domain [0, s] at u and its final 
local minimum at v (see figure |l|). We begin by choosing 



Then setting q = in (|5.2| ), we have 

rf ) =# + ^'(4 J ' ) ), (5-5) 
(i) (i) 

so there is a unique x with v < x < s such that if i/^o = x then =0, and this x 
is independent of j. As increases, x moves towards s independently of j. Then we can 
choose 6 sufficiently large that 

\V'(x)\ < V'(u). (5.6) 

This is the final constraint on 6, so our critical value 6 C will indeed be independent of j. 

(i) (i) 
Now if we increase ipQ from x towards s, i/j^ increases to s. Therefore, we can choose 

y such that x < y < s and when ip^ = y then ip^ = u. 



Using (|5.2|), each ijj^ defines a sequence ipi . We next define the set 

A® = |4 J) G [x,y) :3q,0<q< 2j,^y'(^' ) ) > o[ . (5.7) 

fc=0 



We see from ( |5.2| ) that if ip^ leads to a sequence ip^ which is not monotonically decreasing, 
then V>o J) 

is in A^> . In addition there may be a monotonically decreasing sequence for which 

4 j) is in A®, if E 2 k=o V '(4 j) ) > 0- 

We note that x ^ A^\ since the sequence which starts at x jumps first to 0, and then 
gets locked into the negative region so continues to decrease. On the other hand, using 
J5.6p , we see that y — e £ A^ as long as e is sufficiently small. So A^ is not empty. Since 
V' is continuous, we also note that A^' is open. 

Following [14], we choose 

4 J) =inf (5.8) 
Then ^ £ A® since is open, and so 



^n^ j) )<0, ? = 0,..-,2j. (5.9) 



fc=0 

But since is the infimum of there must by continuity be at least one q such that 



and we choose q to be the smallest such q 



EVrf) = (5-10) 
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Let us suppose q < 2j. Then since X^fe=o V'wk ) < ^, we have certainly V'(ip~ ) > 0. 

9+1 ^9 



Using 0, = and therefore EjK^'W^) > °> which contradicts Q. We 



deduce that Efcio = 0, and so we have completed the construction of a solution 

to the equations (|5,2|), fl5,3|) for sufficiently large 9. This solution has the property that the 



sequence ip k is monotonically decreasing. 

We can extend the proof given above (in the same manner as in [15]) to construct 
further soliton solutions. This time, we fix an integer N with < N < 2j and require 

V'[v) 

Then there is a unique x& such that if i()q = then 

ij$ > # > ■ ■ ■> V#Li > v (5-12) 

and V? = 0. 

Only in the case of N = 1 which we covered before is x^' independent of j. However, 
we note that x^> > x^ +1 \ and define linij^oo x^' = x^°°^ £ (v,s). If we take the limit of 
equation ( |5.2| ) as j — ► oo, we see that as we increase to infinity, a^ 00 ) increases to s. This 
means that we can make each x^' arbitrarily close to s by increasing 9 independently of j. 

We next take 9 sufficiently large that 



N 



y'(£ (oo) ) <V'(u). (5.13) 



This is the final constraint on 9 and we note that it depends on N, but not j. We choose 
y(j) > %(J) such that if ip^ = y^ then V'jv+i = u - Replacing x and y with x^ and in 

the definition ( |5.7D of the proof now proceeds as before. 

— (i) 
We have constructed for each N = 1, 2, • • ■ , 2j a soliton with the property that the ip q 

form a monotonically decreasing sequence. The first N eigenvalues lie in the interval (v, s) 

and the others lie in (0, u). In particular, we note that each tpq satisfies V"(i/j^) > 0. In 

the large 9 limit, the first N eigenvalues tend to s and the others to zero. 



6. Stability 

In this section, we investigate the stability of diagonal solitons. We will show that a 
necessary condition for stability is that the eigenvalues form a monotonic sequence. Next, 
we specialise to values of fx greater than some critical value \x\ depending on the potential 
V and on j. Then another necessary condition for stability is that the eigenvalues lie in 
the regions close to the troughs of V in which V" > 0. Together, these two conditions are 
also sufficient for stability, so that the solitons we constructed in section |5| are stable for 
sufficiently large [i. Finally, we show that the only deformations of a stable soliton which 
leave its energy invariant are rotations. 
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The tools for analysing stability have been set up and discussed in detail in [15] in the 
context of the Moyal plane (and see also [24] in which stability under radial perturbations 
was considered). The stability functional is defined 



2 



(6.1) 



where (j) is a soliton solution. If E is nonnegative then the soliton (p is said to be stable. 
We will find that E always has zero eigenvalues because of the rotation symmetry of the 
action. 

Substituting for the action, we can rewrite the stability functional as 



d 2 

E(w) = J(w) + fi-^ Tr V{4> + euj) 



(6.2) 

e=0 



where J(oS) = Tr [Jj,u;][u;, Jj\ and /i = R 2 /2. Using standard non-degenerate perturbation 
theory (see [15]), this can be written 

EM = J(u>) + 2/x V |(n|u;|m)| 2 T//( ^ ) ~ V ' { ^ n) + /i V |(n|o;|n)| V'(^) (6.3) 

as long as m ^ 0^ for all /n, n. If there are degenerate eigenvalues <pm — (ftni we can use 
the fact that the expression (|6.2j ) is continuous in (f>, and obtain the correct formula by 
taking the limit of ( |6.3[) as the eigenvalues approach each other. In this way, we see that 
we must make the substitution 

V'(^ m ) - V'jM 

7 7 ► V (0 m ) (6.4) 

4>m ~ <Pn 

in equation ( |6 . 3| ) . 

Examining the equation of motion ( |2.6| ), we see that the eigenvalues must lie within 
0(/i _1 ) of the zeros of V . If any </> m is close to r, then since V"(r) < 0, the last term in 
Q6.3| ) can be made large and negative, so the soliton is unstable. So choosing fi sufficiently 
large, our first necessary condition for stability is that each <p m lies in the region [0, u] or 
[v,s] (see figure [l]). This condition simply says that the eigenvalues must all lie in the 
regions close to the troughs of V in which V" > 0. 

We can write 

J(oj) = \yj a n -i{n\u:\m) — ^/a m _i(n — l|w|m — 1)| 2 + (n — m) 2 |(n|(j|m)| 2 (6-5) 

m,n 

and note that E is a quadratic form in the matrix elements of uj, and these matrix elements 
are only coupled along diagonals, which we label by an integer k. We dealt with the case 
k = in the previous paragraph, and so the stability problem can be reduced to studying 
the reduced functional 

E fc (w) = ^2 \^a n -i{n\uj\m) - yja m -i{n - l\u\m - 1)| 2 

\n—m\=k 
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+ (n — m) 2 \ (n\uj\m)\ 2 

\n—m\=k 

+2fi > nuh T 7 ( 6 - 6 ) 

^ , <Pm-<Pn 
n— m=fe 

for each A; > 0. 

We fix A; and define x n = (n + /c|o;|n) for n = —j, ■ ■ ■ , j — k, so that x* = (n|<jj|n + A;). 
Then we can write 

£ fc (a;) = 2 ^ gwx m x* (6.7) 



m.n 



where the only non-zero elements of the quadratic form q are 

Qnn = a n + a n+fc _i + fc(A: - 1) + -y n (6.8) 

gn,n-l = -yJCtn-lOln+k-l (6.9) 

and where 

In = Ai 7 7 > <?n 7^ <Pn+fc (6.10) 

<Pn+fc — <?n 

7 n = ^V"{<p n ) , 0„ = n+ fc. (6.11) 

To obtain (|6.8|), we used the identity 

7j(an+fc-i + Oi n + a n _i + a n+fc ) = a n + a n+k ~i - k. (6-12) 

We will follow [15] to use elementary row and column operations to find a new diagonal 
quadratic form, C, which has the same numbers of positive, negative and zero eigenvalues 
as q. Specifically, we define the diagonal elements of C inductively by 

C-i = Q-j-j (6-13) 
C„, = q n ,n ~ , (n = -j + 1,- ■■ ,j - k). (6.14) 

We are now ready to prove another necessary condition for stability. 

Theorem 2 A necessary condition for stability is that the eigenvalues 4> n form a mono- 
tonic sequence. 

We note first that we do not need to assume /j, is large. 

We examine Xfc(u;) with k = 1. We can use ( |2.9[) to eliminate V'(4> n ) and V'((j) n +i) from 
7 n , and find 

4>n+2 — 4>n+l . 4>n ~ 4>n-l , a -. r\ 

qnn = a n+ i— hOn-i- — , (6.15) 

Pn+l ~ <Pn 4>n+l ~ <Pn 

as long as <f> n ^ <t>n+i- Let us assume first that there are no such pairs of consecutive 
degenerate eigenvalues. Then it is very easy to check by induction that 

„ 4>n+2 - 4>n+l ,„ 1c v 

C n = a n+ i— — (6.16) 

4>n+l ~ <Pn 
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for all n = —j, ■ ■ ■ , j — 1. If the (ft n are not monotonic, at least one of the C n must be 
negative. 

For the degenerate case, assume that the first pair of consecutive degenerate eigenvalues 
is (f) p = (ftp+i- Then we obtain 

C n = a n+ i— — , n<p-l, (6.17) 

9n+l ~ <Pn 

and in particular, C p _i = so that this inductive reduction fails at n = p — 1. However, 
we can use the (p — l)-th row and column to further eliminate and bring the matrix into 
block diagonal form, and then C p and C p+ i are given by the eigenvalues of 



-^/a p _ia p _ 2 
-^ap-xoip-2 1a p + fJ-V"((p p ) 



(6.18) 

one of which is always negative. This completes the proof of Theorem 2. 



In equation ( |6.16|) , we note that Cj-i is zero (since ay = 0). This zero eigenvalue in 
the stability functional comes about because we can rotate the soliton, and obtain a new 
solution with the same energy. In general, if there is a family of soliton solutions 
depending smoothly on v, with <p(0) = (ft, then necessarily there exists an w ^ with 
S(w) = 0. To see this, take u = d(j)/dv\ v= Q. We can always choose a parameterization such 
that Then, expanding (j) in terms of some basis for the matrices, 4> = </>*T*, we have 



2 



o = i^[*(,)] 



i d 2 s 

u=0 ~ 2 dtfdt 



dS d 



•2 Mi 



(6.19) 

v=0 



The first term is the stability functional, whilst the second term is zero by the equation of 
motion, so we obtain T,(uj) = 0. 

Since the rotations are a symmetry of the action, they will lead to zero eigenvalues 
of the stability functional. The rotations are a three parameter group, but one of these 
(rotation about the z-axis) leaves a diagonal soliton invariant. This means we would expect 
to find two zero eigenvalues in the stability functional corresponding to the rotations. The 
quadratic form q acts on a complex space, and so the zero eigenvalue Cj_i corresponds to 
two zero eigenvalues on a real space. Our task now is to show that this zero eigenvalue 
corresponds precisely to the rotations. 

We do this by showing that any hermitian matrix X £ Mj can always be rotated into a 
canonical form. 1 Then we can restrict our space to matrices of this form, and check that the 
zero eigenvalues disappear. Specifically, given any hermitian matrix X pq , (—j < p,q < j), 
we can always use a rotation to set Xjj—i = = 0, and we shall restrict to this 

rotation fixed space of matrices. 

To see that this is possible, it is helpful to use the coherent state representation for 
the algebra. The formulae that we use can be found in [36]. A general rotation generated 



1 We recall this notation from section g. The fuzzy sphere is the algebra Mj of (2j + 1) x (2j + 1) complex 
matrices. 
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by J\ and J 2 can be re-written in normal form 



exp(i#i Ji + i9 2 J 2 ) = exp(C J+) exp^Jg) exp(-CJ-) = D(C) (6.20) 

where, defining 6> = #1 + z#2> we have £ = tan(|0|/2) and 77 = log(l + |C| 2 )- It can 

also be written in antinormal form 

D(C) = exp(-CJ-) exp(-r ? J 3 ) exp(CJ+). (6.21) 

The coherent state is defined 

\z) = D(z)\j,-j), (6.22) 

and is automatically normalised (z\z) = 1. 
The operators D satisfy 

D((i)D(( 2 ) = £>(Ca) exp(id> J 3 ) (6.23) 

where £3 and $ are respectively complex and real numbers depending on (j and £2- I n 
particular, D(£)D(— Q = 1 as can be checked by multiplying the normal and antinormal 
forms. 

Given a hermitian matrix X E Mj, we define the covariant symbol [37] 

X(z,z) = (z\X\z), (6.24) 
which, in terms of matrix elements in the standard basis, is 



*^>=(TTm^ £.*™^" + VUJUJ- (6 ' 25) 

V 11/ m,n=—j V v " y \ •> s 

The covariant symbol is one way of mapping the fuzzy sphere algebra Mj to an algebra 
of functions on the sphere (with noncommutative *-product). Here the sphere is given in 
stereographic coordinates, and we note that if the matrix X is diagonal, X(z,z) depends 
only on \z\, so that it is indeed rotationally invariant about the north-south axis. 

The function X(z,z) must have at least one stationary point, zq, in \z\ < 00. One can 
quickly check using (HH) that applying a rotation X -> X = D\C,)XD(Q) leads to a new 
covariant symbol X(z,~z) = X(z',~z') for some z' depending on z and £. In particular, we 
can apply the rotation D{— Zq) to get a covariant symbol with a stationary point at z = 0. 
Then differentiating ( 6.25|) at z = gives 



Xjj-i = Xj-tj = (6.26) 

which is our desired rotation fixing condition. 

We can now deal with the zero eigenvalue Cj-i- We restrict to the rotation fixed space 
of matrices X, with Xjj-i = Xj_ij = 0. In this case, the above analysis works identically 
except that Cj—i no longer appears. This shows that the zero eigenvalue corresponds 
precisely to the rotations. In particular, when the eigenvalues form a monotonic sequence, 
Si is positive definite on the rotation fixed space. 
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We now go on and consider for k > 1. 



Theorem 3 Suppose <fr is a diagonal soliton, and the eigenvalues 4> m form a strictly mono- 
tonic sequence, and suppose /i is sufficiently large. Assume further that every eigenvalue 
4> m has V"{<pm) > 0. Then the quadratic forms are positive definite on the rotation 
fixed space. 

We have already covered the cases k = and k = 1, so we restrict now to k > 1. Since \i 
is large, the condition V"(<f> m ) > means that the eigenvalues can lie close to and s, but 
not r. We choose, without loss of generality, the sequence of eigenvalues to be decreasing 
(rather than increasing). Then there exists N such that 

<!>-,, ■ • • , <f>- j+N ~l = s - 0{^ 1 ) (6.27) 
<f>- j+N ,--.,<t>j = 0(ii- 1 ) (6.28) 

that is the first N eigenvalues are close to s, and the remaining are close to 0. We can use 
(2.10) to compute the 4> m to first order in fi~ x . We see immediately that if 4> m = s—0(fi~ p ), 



then (f)m+i = s — 0(n p+1 ) for m < —j + N — l. Then we can write 

-j+N-l 

-S + 0(lT l ) = 4>-j+N - (j>-j+N-l = Y] V\(j) k ) 

<*-j+N-l 

11 -V'^+n^ + O^- 1 ) 



giving 



and 



Similarly, 



and 



at-j+N-i 

^ F"( S )(^- +J v-i -s) + Oifi- 1 ) (6.29) 

j+N-l 



-j+N-i = s - '""^"V 1 + 0(^ 2 ) (6.30) 
cp- j+N - q = s - 0{ix- q ) , q>2. (6.31) 



^+N= Sa v i^- 1 ^ 1 + 0(^) (6.32) 



-j+N-i+ q = 0(fi- q ) , q>2. (6.33) 



We can use these results to evaluate the 7 m to first order, and find bounds on the C m . 
First, when m + k < — j + N — 1, we have j m = /j,V"(s) + 0(1) giving 

C m = iiV"(s) + 0(l), m<-j + N-k-l. (6.34) 

Next, if m = — j + N — k, we find 7-j+jv-fc = —ct-j+N-i + 0(m -1 )) giving 

C^+N-k = ot- j+N - k + k{k - 1) + Oifi- 1 ). (6.35) 
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Then for —j + N — k < m < —j + N — 1, we have — CO -1 ), and we can check 
inductively that 

C m > a m + k{k-l) + 0{^ 1 ), -j + N-k<m<-j + N-l. (6.36) 
When m = — j + N — 1, we have 7_.7-4.jv-1 = — a-j+jv-i + 0{[i~ l ) giving 

C-j+N-i > k(k - 1) + ©(/i" 1 ). (6.37) 
Finally, if m > — j + N - 1, we have 7 m = nV"(0) + 0(1), and so 

C m = fiV"(0) + 0(1) , m>-j + iV-l. (6.38) 
This completes the proof of theorem 3, since every C n is positive at sufficiently large fi. 



We have shown that a soliton which satisfies the assumptions of theorem 3 has positive 
definite stability functional on the rotation fixed space, and so is stable. In particular, the 
solitons which we constructed in theorem 1 are stable when \i is sufficiently large. These 
solitons cannot be deformed without increasing their energy, because this would require 
the stability functional to have zero eigenvalues after fixing the rotations. 



7. Unstable multi-lump solutions 

In this section we consider solutions of the field equations that describe separated lumps 
on the fuzzy sphere. As mentioned in the Introduction, a precise mathematical definition 
of a multi-soliton is lacking at finite non-commutativity but one can nevertheless look 
for solutions where the scalar field has nonvanishing values in well separated regions on 
the sphere. Perturbative calculations at large non-commutativity indicate that such lumps 
attract each other [12,13] and therefore one would not expect to have solutions with multiple 
lumps at arbritrary locations. On symmetry grounds, however, one might anticipate that 
there could exist solutions with lumps in an unstable equilibrium. The simplest such 
solution would have two identical lumps, one at each pole of the sphere. Such a soliton 
would be rotationally invariant about the z-axis but also invariant under reflection about 
the equatorial plane. It would thus be a diagonal solution with the eigenvalues satisfying 



(7.1) 



If j is a half integer, the equations ( 2.1C| ), ( |2.11 ) reduce consistently under this assumption 
to 

R 2 m 3 
b m +i-4> m = - y^V'i^i), m = -j, (7.2) 

*=-3 



with the constraint 
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To see this, we use the property a_ m = a m -l to check that ( 2 . 1 0|) also holds for m > ^. 
When j is an integer, the equations become 



R 2 m 

1 ~ 1 = ~ ^ ^'Ofo) » m = -3, ■ ■ ■ > - 1 ( 7 - 4 ) 



2o! m . 



with the constraint 



-l 

2 F'(^) + V'(4>o) = 0. (7.5) 
i=-j 

We can suitably adjust the set (equation (|5.7|) ) to the altered difference equation and 
constraint. Then applying the argument of section ||] shows immediately that these soliton 
solutions do exist for sufficiently large 9. We see from the stability analysis in section || 
that such solitons are unstable since the eigenvalues do not form a monotonic sequence. 

We can construct a bifurcation diagram for these solutions which must be a subset of 
the bifurcation diagram for the full problem. Then we see that they correspond precisely 
with certain trajectories in the full bifurcation diagram. In particular, the diagonal eigen- 
states of [Jj, [Jj, • ] ] corresponding to eigenvalue k(k + 1) with k even have the property 
<ftm = <ft-m, so every second bifurcation from x = r corresponds to an unstable solution of 
this type. 

By symmetry, one would expect to be able to go further, and exhibit unstable solitons 
with n lumps evenly spaced around a great circle. At fixed j, this n can be at most 2j, 
as the fuzziness prevents resolution of more lumps. Here, we briefly report results for a 
specific example. We examine j = § and seek a 3-lump solution. We begin with a general 
diagonal matrix P, and a rotation through 2tt/2> about the x-axis, 1Z = exp(27ri,J x /3), and 
define a new matrix Q 

Q = P + KPK- 1 + U- l PR (7.6) 

which is invariant under rotation through 2tt/2> about the x-axis. In this case, the situation 
is simplified because all such matrices Q, and kinetic terms [Jj, [Jj,Q]] are simultaneously 
diagonalisable, and two of the eigenvalues are always identical. One can then solve the 
equations of motion for the eigenvalues, and find several solutions when fi is chosen large 
enough. 

In the following section we take a limit of large R and j to blow up a small region 
close to the south pole into the Moyal plane. The multi-lump solutions constructed in this 
section will not survive this limit since any lump not sitting at the south pole will be moved 
to infinity. However, by symmetry one would expect the existence of an unstable solution 
consisting of an infinite chain of lumps finitely separated along a straight line in the Moyal 
plane, and its generalisation to a two-dimensional array of lumps. 2 Such solitons would of 
course have infinite energy. 

8. The planar limit 

An additional motivation for studying the fuzzy sphere is that, by blowing up a small 
region close to, say, the south pole, there is a weak convergence of the algebra to that of 
2 In [16] such periodic solutions are constructed in the limit of large 8. 
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the Moyal plane [19,34]. This means that the fuzzy sphere may be used to cut off the 
infinite dimensional algebra and provide an infra-red regularisation for field theory on the 
Moyal plane, but care is required when taking advantage of this limit because of the global 
differences between the sphere and the plane. In this section we show that the solitons 
constructed in section || do converge in a weak sense to solitons on the Moyal plane. The 
relevant limit is to send j and R to oo, keeping the ratio 

8 = — (8.1) 

fixed. In theorem 1 of section @, we found diagonal solitons for 9 large enough. For each 
< N < 2j, we constructed a soliton ipj: with the first TV" eigenvalues lying in the interval 

(v,s) and with the others in (0, u) (see figure pi). The constraint on 9 is independent of j. 

(i) 

We will show in theorem 4 below that the solutions tpY have the property that there exist 
tpk with 



(a) 

W k Wk as j oo 

at fixed k. Assuming this result for the moment, we can take the limits of equation (| 
at fixed q and discover 

9 q 

^ +1 -^=— -VV'W, q = 0,l,--- (8.3) 



2(q + l 



k=0 



which is precisely the equation for a diagonal soliton on the Moyal plane with noncom- 
mutativity parameter 8 [1, 14]. We will in addition see that ipk — > as k — > oo which is 
necessary for the Moyal plane soliton to have finite energy. 3 

Theorem 4 Send j — ► oo and R — > oo keeping 9 fixed and sufficiently large. In this limit, 
the matrix elements of the fuzzy sphere solitons constructed in section || converge to the 
matrix elements of a soliton on the Moyal plane. 

We begin by showing that ip^ = inf is a decreasing sequence and so has a limit. For 
each j, we define the interval i, = [x^\y^'], where x^\y^' are the points defined in the 
proof of theorem 1, so that A^' C Ij. We recall that x^ +l ^ < x^> and y^ +l ^ < y^>. 
Therefore, if Ij D Ij+i = 0, we have inf A^ +1 ^ < inf A^ and we are done. 

If Ij n Ij + \ / 0, then pick a £ Ij Pi Ij + \ and suppose a £ A^ +1 \ If we set rp^ = 
ipQ +1 ^ = a, then a ^ A^ +l ^ means that 



k=0 

We claim that 



£n^r J )<0, Vq- (8.4) 


< , Vfc. (8.5) 



3 We note here that the stability analysis for the fuzzy sphere does not carry over to the Moyal plane by 
simply taking this limit. The condition given on R for stability does not match the condition on 9 for the 
planar limit, and the more technically complicated analysis given in [15] is therefore still required for the 
Moyal plane. 
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By assumption, this is true for k = 0. We assume that (|8.q ) is true for k < n and argue by 

and V?" 



induction. Since, for fixed k, and ipft both lie in a domain in which V' is increasing, 



we have 

V'ty®) < V'{^ +1) ) (8.6) 

for k < n. Then 



+ W 2(j+l)i fc=0 



- 2(n + l)(l-^ Ty ) 



k=0 



it 



and, by the induction hypothesis, this implies that 



^"^■ W + 2(„ + l) (1 - § ) g^ 



establishing Next, using (^J) and (|8J3|) , we have 



9 



^% W )<0 Vg. (8.9) 

fc=0 

which means a £ 

So far, we have shown that if a ^ .AW+l) then a ^ j4y°) and this implies 

n c ij n c A^ +l \ (8.10) 

Taking the infimum, 

inf(/ i+1 n A&) = inf > inf A {j+l \ (8.11) 

so is a bounded, decreasing sequence, and so has a limit, ip^ — > ipo, as 

j oo. 

The next stage in proving theorem 4 is to show that, for each fixed k, there is a V'fc 
such that Vfc ~~ * Vk as J ~ *• 00 • This follows easily by induction. Assume it is true for 
k < n and use (15.21) to note 



2(n + 1)(1 - £) 

n 

-^ + ^— ryE^) ( 8 - 13 ) 

V ; fc=0 

as j — * oo by continuity of V'. 



-21- 



Finally, it remains to check the finite energy condition tp^ —* as k —* oo. We first 

note that, since any soliton on the fuzzy sphere has its spectrum in the interval [0, s], we 

(?) (i) (i) 

have > for all k,j, and this implies ^ > for all k. Similarly, since — V'Jt > 

we also have V'fc+i — "4>k- Then V'fc+i is a decreasing sequence which is bounded below and 

so there exists > with ^ — ► ^ as k — > oo. 

(7) 

By our original construction, we also have < u for all j and k > N, and so, taking 
the limit, ipk < u for k > N. Then only ipo, ■ ■ ■ , tpN he outside the region in which V' is 
increasing, so V'(tpk) > V'(t/)) for k > N + 1. Substituting into fl8.3Q gives 



and we deduce that ip = 0. So ^fc — > as A; — > 00, and this completes the proof of theorem 4. 
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